On the ACDM Universe in f{G) gravity 



X 



Ratbay Myrzakulov 1 , Diego Saez-Gomez 2 and Anca Tureanu 3 
1 Dept. Gen. Theor. Phys., Eurasian National University, Astana, 010008, Kazakhstan 
2 Institut de Ciencies de I'Espai ICE/CSIC-IEEC, Campus UAB, 
Facultat de Ciencies, TorreC5-Parell-2a pl,E-08193 Bellaterra (Barcelona) Spain and 
3 Helsinki Institute of Physics, P.O. Box 64, FI-00014 Helsinki, Finland 

In the context of the so-called Gauss-Bonnet gravity, where the gravitational action includes 
function of the Gauss-Bonnet invariant, we study cosmological solutions, especially the well-known 
ACDM model. It is shown that the dark energy contribution and even the inflationary epoch can be 
explained in the frame of this kind of theories with no need of any other kind of component. Other 
cosmological solutions are constructed and the rich properties that this kind of theories provide are 
explored. 
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I. INTRODUCTION 
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It is widely accepted by the scientific community that our Universe is currently in an accelerated phase. General 
QO ■ Relativity in its standard form can not explain the accelerated expansion without extra terms or components, which 
have been gathered under the name of dark energy. Since the discovery of the acceleration, a large number of possible 
mechanisms have been proposed to explain the origin of the dark energy, from the cosmological constant, scalar fields, 
■ to modifications of general relativity as well as other alternatives (for reviews on unified inflation-dark energy modified 
q-i gravities, see Refs. [Hfl], an d for comparison with observational data see Refs. [Ij-O)- The most popular idea is 
represented by the cosmological constant, whose origin may be explained by means of the vacuum energy density, 
although its value does not match the one predicted by quantum field theories. Even in such a case, the Equation of 
State (EoS) parameter is a constant equal to —1, while some recent observations suggest that the EoS parameter is a 
dynamical variable and it could have crossed the phantom barrier, with the EoS parameter being less than — 1. Such 
a dynamical behavior of the EoS could be well explained by scalar fields with quintessence or phantom behavior (see 
^vq ' Ref. [5]). Another alternative is the modification of the gravity law. Several ways have been suggested to perform 
such a modification of General Relativity by means of the modification of the Hilbcrt-Einstein action, which obviously 
0^ ' modifies the field equations. In this framework, the so-called F(R) gravity has been explored, which introduces a more 
complex function of the Ricci scalar in the action and explains well the cosmic history via cosmological reconstruction 
q{ (see Refs. @-[Hj]). 

In this work we study the so-called Gauss-Bonnet gravity, where the gravitational action includes functions of the 
Gauss-Bonnet invariant. This kind of theories have been investigated and may reproduce the cosmic history (see Refs. 

[36[). Here, we show that the ACDM model can be well explained with no need of a cosmological constant but 
with the inclusion of terms depending on the Gauss-Bonnet invariant in the action. Even more, it is shown that the 
extra terms in the action coming from the modification of gravity could behave relaxing the vacuum energy density, 
represented by a cosmological constant, and may resolve the so-called cosmological constant problem. To study and 
reconstruct the theory that reproduces such a model as well as other kind of solutions studied, we shall use a method 
proposed in Ref. § for f(R) gravity and implemented for /(G) gravity in Ref. [36j], where the FLRW equations 
are written as functions of the so-called number of e-foldings instead of the cosmic time. The possible phantom 
epoch produced by this kind of theories is also explored, as well as other interesting cosmological solutions, where the 
inclusion of other contributions as perfect fluids with inhomogeneous EoS are studied. 

II. [R + f(G)] GRAVITY 

We consider the following action, which describes General Relativity plus a function of the Gauss-Bonnet term (see 
Refs. and 0): 



S = Id Xy/—g 



-l 1 R + f(G) + L r , 



(i) 



where k 2 — 8itGn, Gm being the Newton constant, and the Gauss-Bonnet invariant is defined as usual: 

G = R 2 - AR^R^ + R^ x *R^ Xa ■ (2) 
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By varying the action over <7 M „, the following field equations are obtained: 

= ^(-^ + \^ VR ) + T ^ + l^fiG) - 2f G RR^ + AfaRpR vp 
-2f G RW° r R v paT - Af G R^R pa + 2(V»V v f G )R - 2g^(\7 2 f G )R - 4(V P V/ G )i?^ 

- 4(V p V' y / G )i? w + 4(V 2 /g)^ + 4 5 ^(V p V ff / G )i?" CT - 4(V P V CT / G )i? ww , (3) 

where we made the notations f G — f'(G) and / GG = f"(G). We shall assume throughout the paper a spatially-flat 
FLRW universe, whose metric is given by 

3 

ds 2 = -dt 2 + a{tf^{dx 1 ) 2 , (4) 
i=i 

where a(t) is the scale factor at cosmological time t. For the metric (|4]), the field equations give the FLRW equations, 
with the form 

= -^H 2 + Gf G - /(G) - 24GH 3 f GG + p m , 

= 8H 2 f G + 16H(H + H 2 )f G + \{2H + 3H 2 ) + f- Gf G +p m . (5) 

The Hubble rate H is here defined by H = a/ a, while the matter energy density p m satisfies the standard continuity 
equation: 

p m + 3H{1 + w)p m = , (6) 

while the Gauss-Bonnet invariant G and the Ricci scalar R can be defined as functions of the Hubble parameter as 

G = 24(HH 2 + H 4 ), R = 6(H + 2H 2 ). (7) 

Let us now rewrite Eq. ([5]) by using a new variable, N = In ^ = — ln(f + z), i.e. the number of e-foldings, instead of 
the cosmological time i, where z is the redshift (this method has been implemented in Ref. Q for f(R) gravity). The 
following expressions are then easily obtained 

n-nr N H - N - dN d - H d ^ - H 2 * + H H> d H> - dH (to 

a-a e , H-N-—, Jf -H—, —-H— + HH-, H- — . (8) 

Eq. ([5]) can thus be expressed as follows 

= -^H 2 + 24H 3 (H' + H)f G - f - 576H 6 iHH" + 3H' 2 + iHH') f GG + p m , (9) 

where G and R are now 

G = 24(H 3 H' + H 4 ) , G = 2A(H 4 H" + 3H 3 H' 2 + AH 4 H') , R = 6{HH' + 2H 2 ) . (10) 
By introducing a new function X (XS X = H 2 we have 

H' = -x-^x' , H" = --x- 3 ' 2 x 12 + -x- l ' 2 x" . (11) 
2 ' 4 2 v 1 



Hence, Eq. © takes the form 

fGG + Pm , (12) 



= -4rx + 12x(x' + 2x)f G - f - 2A 2 x 



^x 2 x" + ^xx' 2 + 2x 2 x' 



where we have used the expressions 

G = Ylxx' + 24x 2 , G = 12x-^ 2 [x 2 x" + xx' 2 +4x 2 x'} and R = 3x' + 12x. (13) 

Then, by using the above reconstruction method, any cosmological solution can be achieved, by introducing the given 
Hubble parameter in the FRW equations, which leads to the corresponding Gauss-Bonnet action. 
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III. RECONSTRUCTING ACDM MODEL IN R + f(G) GRAVITY 

We are now interested to reconstruct ACDM solution in R + /(G) gravity for different kind of matter contributions. 
It is well known that such a solution can be achieved in GR by introducing a cosmological constant (cc) in the action. 
Nevertheless, we show that in GB gravity there is no need of a cc. The cosmological models coming from the different 
versions of modified GB gravity considered will be carefully investigated with the help of several particular examples 
where calculations can be carried out explicitly. 

In this paper, we restrict ourself to explore some "classical" modified gravity models for the ACDM case, as well 
as other interesting and important cosmological solutions. For the ACDM model, the Hubble rate is given by 

» 2 = t + 0' (»> 

where po is the matter density (which consists of barionic matter and cold dark matter) and A is the cosmological 
constant. In the rest of this section we put K = 1. 

For the ACDM model, described by the Hubble parameter (fT^fj) , we can write the derivatives of the scale factor as 
well as the Hubble parameter in the following useful way: 



Aa 2 p Q .. 2Aa 3 - p 
3 3a ' 6a 2 



2a 3 2 V 3 J ' 2a 3 V 3 3a 3 2 V 3 

Using these formulas we get 



Then, the following relation between R and G holds: 

G = -^(i? 2 -9Ai?+18A 2 ) . (17) 
Let us recall that x = H 2 . Then, some of the above formulas take the form: 

H=^(A-3x) H = ^(3x-A)y/x~ , R = Z(A + x), 

G = 12x(A — x) , R = 3(A - 3x}y/x , G = 12(3x - A) (2a; - A)y/x . (18) 
Note that the variable x can be expressed in terms of R or G as 



R , 3A ± V9A 2 - 3G . , 

x = - - A or x= , 19 

3 o 

respectively. The above formulas will be useful to reconstruct the ACDM model as well as other cosmological solutions 
in the context of Gauss-Bonnet gravity, as it is shown below. 
We write the first Friedmann equation ([5]) in the form 

= -'iH 2 + 12H 2 {A - H 2 )f G -f- 288H 4 {3H 2 - A)(2H 2 - A)f GG + p m , (20) 

or 

= -3a; + 12x(A - x)f G - / - 288x 2 (3:e - A) (2a: - A) f GG + p m . (21) 
For further algebra more convenient is the following form of this equation 

0={ Pm -3x- /)(A - 2x) + [48x 2 (3x - A) + a; (A - x)]f x + 2Ax 2 {3x - A) (A - 2x)f xx . (22) 
Now we wish to construct some particular exact solutions of this equation. 



A. Case I: p m = 

First of all, let us consider the simple case in absence of matter, p m = 0. Then, the equation (|22l) takes the form 
= -(3a; + /)(A - 2x) + [48a; 2 (3a; - A) + x(A ~ x)]f x + 24x 2 (3a; - A)(A - 2x)f xx . (23) 
We can analyze the cases where the cosmological constant term in the solution (|14p vanishes and where it is non-zero. 
• Let A = 0. Then Eq. ([23]) reads as 

= -2(3a; + /) - x(lUx - l)f x + lUx 3 f xx . (24) 
The general solution of (|24"|) is given by 

f(x) = C 2 x 2 + C*ix(144a; - l)eiw + Vl (x) , (25) 



where 



Here 



Vl (x) = -864a; 



+ xlnx+( — — x ) Ei ( 1, - — I e'"» 



144a; \ 144x / V ' 144a; 



(26) 



/OO 
e- zs s" a ds . (27) 



The two integrals of motion are 



_ 144 
1 x 



.<■(! - 144r)(/, - ) i ( 288a;-2 + ^ ) (/ - n) 



(28) 



J 2 = r— 1 ~ — {f-vi- 144[a;(l - U4x)(f x - v lx )+ 
a;(144a; — l)ei«^ 

288 "~ 2+ 14la^) (/ " Wl)] } ' (29) 

• Let A 7^ 0. Then Eq. (|23p has a complex solution, which has no physical meaning as it gives a complex action. 

Hence, it appears that the ACDM model (fT4|) can not be reproduced by R + f(G) gravity in the absence of matter. 
The only solution found restricts the Hubble rate to give a decelerated Universe. 

B. Case II: p m ^ and A = 

We now explore the case when some kind of matter with a particular EoS is present in the Universe, but with no 
cosmological constant term in the Hubble parameter described in (|T4)) . We explore several examples where different 
kind of matter contributions are considered. 

1. Example 1 

Let us now consider the case when A = and the evolution of the matter density behaves as 

p m = 3H 2 = 3x. (30) 
In this case the modified Friedmann equation ([5]) reads as 

= 2/ + x(144a;-l)/ x -144a; 3 /^- (31) 
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This equation has two integrals of motion 



h = 



h 



1 



144 



s(l - 144a:)/* + ( 288a; - 2 + -j- 



144a; 



/ - 144 



x(l - 144.x)/, + 288a; - 2 + — — / 



X(1AAX — l)ei44a: 

The general solution of the equation pip is given by 

f(x) = C lX 2 + C 2 x(lUx - l) e ™r. 
This function reproduces the solution (fl~4f under the conditions imposed above. 



144a; 



(32) 
(33) 

(34) 



2. Example 2 

Now we consider a more general case, where the energy density is given by, 

p m =u(x) , 

where u(x) is some function of x. In this case the modified Friedmann equation ([5]) reads as 

= 2[3x - u(x) + /] + x(lUx - l)f x - lUx 3 f xx . 

Its general solution is 

f(x) = Cix 2 + C 2 x(l44x - l)e™s + v 2 (x) , 

with 



v 2 {x) = 288x 



1 \ _J_ x 

x - — I e 144 ^ Ji - J 2 



144 



where 



Ji 



3a; — u i 



e 144a = o?x, J2 



/ 



144 



(3a: -u) (144a; - 1) 



Note that the two integrals of motion are 



144 1 

I = [ x 1 - 144a; - U2x + (288a; - 2 + — — )(/ - v 2 )] 

x 144a; 



— {/- v 2 - 144 [a;(l - lUx)(f x - v 2l 



x(14Ax - l)e" 



In fact, we can directly verify that 



{ 288x - 2+ iL ]ti 



/la: — — 



(35) 
(36) 
(37) 

(38) 
(39) 
(40) 



(41) 



(42) 



Then, the solution (|3T[) gives the function of the Gauss-Bonnet invariant that reproduces this model for any kind of 
EoS matter fluid. 



C. Case III: p m and A / 

Let us now explore the most general case for the solution ([M)) in R + f(G) gravity with a non vanishing matter 
fluid with a given EoS parameter. 
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1. Example 1 

We consider 

p m = 3x + p. (43) 

Then Eq. ([5]) takes the form 

= (P ~ /)( A - 2ar) + [48x 2 (3a; - A) + a; (A - x)]f x + 24x 2 (3x - A) (A - 2a;)/x* (44) 
and has the following particular solution: 

/(:r) = 7a; 2 - jAx + /3. (45) 

If A = 0, then f(x) — ^x 2 + /3. Also if 7 = 0, then the solution takes the form / = /3, which corresponds to the 
cosmological constant. Note that if /? = —A then 

p m = ^f = 3z-A, /(x) = 7a; 2 — 7 Ax — A . (46) 

This gives a solution where the cosmological constant is corrected by the contribution from /(G), what may resolve 
the cosmological constant problem. 

2. Example 2 

Now we consider the density of the energy which is given by 

p m = kx 2 + 3x + (3. (47) 
The corresponding first Friedmann equation reads as 

= (kx 2 +p- /)(A - 2x) + [48ir 2 (3ir - A) + x(A - x)]f x + 24:r 2 (3ir - A) (A - 2x)f xx . (48) 
This equation has the following particular solution: 

f(x)=jx 2 +(J^- 7 A^x + p, (49) 

where we must have A = — sj. As well as in the above example, the given function f(G) produces a relaxation on the 
cosmological constant, which can be seen as the possible resolution of the cc problem. 

IV. COSMOLOGY IN R + f(G) GRAVITY WITH THE PRESENCE OF AN INHOMOGENEOUS FLUID 

Let us now consider the theory described by the action (JTJ) in the presence of a perfect fluid, whose EoS is given by 
the general expression: 

p = w(a)p + r)(a) , (50) 

where w(a) and 77(a) are arbitrary functions of the scale factor. This kind of EoS could correspond to a dynamical 
viscous fluid or possibly the effective EoS that accounts for the extra terms in the gravitational action, as curvature 
terms or scalar and vector fields (see Ref. [l3[ ). Then, the FRW equations ((5$ are written now as 

3H 2 = p + p f{G) , 2H + 3H 2 = -{p + Pf(G) ). (51) 

Here the energy and pressure densities Pf(G) an d Pf(G) are properly defined to account for the extra terms introduced 
by the f(G) function in the action (JTJ), and are given by 

P f (G) =Gf G -f- 2AGH 3 f GG , 
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Pf(G) = 8H 2 f G + 16H(H + H 2 )f G + f-Gf a . (52) 

Then, by combining the FRW equations (|5ip and using the EoS (f5U|) . we could write the inhomogeneous term 77(a) 
as a function of the Hubble parameter: 

T 1 (a) = -w(a)p f{G) + Pf{G) +2H + 3H 2 {l + w{a)) . (53) 

Hence, by specifying a cosmic solution H(a), any cosmological history can be reconstructed in this framework, where 
the contributions from /(G) and the inhomogeneous fluid described by (I5U)) drives the Universe evolution. To show 
this, let us consider the example 

3H 2 = H 2 + H ia m , (54) 

where Hq, Hi and m are constants. This solution reproduces a Universe dominated by an effective cosmological 
constant at the current time and which enters a phantom epoch in the future (as some observations suggest). We 
could consider w(a) — and the f{G) function calculated in (f2"5"j). which reproduces an effective cosmological constant. 
Then, the effective fluid ([5l)|) describes dust-matter at the beginning, when the cosmological constant dominates and 
then it drives the Universe into a phantom epoch, which probably finishes at the so-called Big Rip singularity. 



V. COSMOLOGICAL SOLUTIONS IN PURE f(G) GRAVITY 

We have studied so far a theory described by the action ([1]), which is given by the usual Hilbert-Einstein term plus 
a function of the Gauss-Bonnet invariant, that is assumed to become important in the dark energy epoch. In this 
section, we are interested to investigated some important cosmic solutions in the frame of a theory described only by 
the Gauss-Bonnet invariant, and whose action is given by 

S = J d 4 x^[f(G) + L m ] . (55) 

In this case, the FRW equations are: 

= Gfo — f — 24GH 3 f GG + Pm 

= 8H 2 f G + 16H(H + H 2 )f G + f-Gf G + Pm . (56) 

We are interested to explore some important solutions from the cosmological point of view, as de Sitter and power 
law expansions. 



A. De Sitter solutions 



De Sitter solutions are described by an exponential expansion of the Universe, where the Hubble parameter and 
the scale factor are given by 

H(t) = H Q ^ a(t) = e Hot , (57) 

where Hq is a constant. This kind of solutions are very important, as the observations suggest that the expansion 
of our Universe behaves approximately as de Sitter. It has been shown in Ref. [12J that de Sitter points are critical 
points in f(R) gravity. It is straightforward to see that this is also the case in f(G) gravity. We can explore the de 
Sitter points admitted by a general /(G) by introducing the solution (|57|) in the first FRW equation given in (|56l) . 
which yields 

= Gof G (G ) - /(Go) . (58) 

Here, Go = 24ffg and we have ignored the contribution of matter. Then, we have reduced the differential equation 
(|56| to an algebraic equation that can be resolved by specifying a function /(G). The de Sitter points are given by the 
positive roots of this equation, which could explain not just the late-time accelerated epoch but also the inflationary 
epoch. The stability of these solutions has to be studied in order to achieve a grateful exit in the case of inflation, 
and future predictions for the current cosmic acceleration. 
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B. Power law solutions 



We are now interested to explore power law solutions for a theory described by the action (|55j) . This kind of 
solutions are very important during the cosmic history as the matter/radiation epochs are described by power law 
expansions, as well as the possible phantom epoch, which can be seen as a special type of these solutions. Let us start 
by studying a Hubble parameter given by 

H(t) = ~ -> a(t) ~ t a , (59) 

where we take a > 1. Then, by introducing the solution (JSHJ) into the first FRW equation (|56")) . it yields the differential 
equation 

= -f(G) + Gf G + ^-/ GG , (60) 
a — 1 

where we have neglected any contribution of matter for simplicity. The equation (|60p is a type of Euler equation, 
whose solution is 

f{G) = dG + C 2 G^ . (61) 

Thus, we have shown that power-law solutions of the type (I59p correspond to actions with powers on the Gauss-Bonnet 
invariant, in a similar way as in f(R) gravity, where power-law solutions correspond to an action with powers on the 
scalar curvature, R (see Ref. [!|). 

Let us now explore another kind of power-law solutions, where the Universe enters a phantom phase and ends in a 
Big Rip singularity. This general class of Hubble parameters may be written as 

H(t) = j?- t , (62) 

where t s is the so-called Rip time, i.e. the time when the future singularity will take place. By inserting the solution 
(l6"2"j) into the first FRW equation (1551) , the equation yields 

4a 2 G 2 

= -f(G) + Gf G (G) - — — , (63) 

1 + a 

which is also a Euler equation, whose solution is given by, 

/(G) = C 1 G + C 2 G 1 ^ . (64) 

Thus, we have showed that power law solution of the type radiation/matter dominated epochs on one side and 
phantom epochs on the other, are well reproduced in pure /(G) gravity, in a similar way as it in f(R) gravity. 



VI. CONCLUSIONS 



We have explored in this paper several cosmological solutions in the frame of Gauss-Bonnet gravity, considering 
specially the case of an action composed of the Hilbert-Einstein action plus a function on the Gauss-Bonnet invariant. 
Also pure /(G) gravity has been considered, as well as the possibility of the implication of inhomogeneous terms in 
the EoS of a perfect fluid, which could contribute together with modified gravity to the late-time acceleration. We 
have shown that the ACDM model can well be explained in this kind of theories, which may give an explanation to 
the cosmological constant problem as the modified gravity terms may act relaxing the vacuum energy density. Other 
kinds of solutions in /(G) gravity have been reconstructed. It has been shown that /(G) gravity could explain the 
dark energy epoch whatever the nature of its EoS, of type quintessence or phantom, and even the inflationary phase. 
More complex cosmological solutions would require numerical analysis, but our analysis of a few simple cases has 
already shown that /(G) gravity accounts for the accelerated epochs and may contribute during the radiation/matter 
dominated eras, and it may explain also the dark matter contributions to the cosmological evolution, what will be 
explored in future works. This kind of modified gravity models which reproduce dark energy and inflation, can be 
modeled as an inhomogeneous fluid with a dynamical equation of state, what would be distinguished from other 
models with a static EoS. Even as perturbations in modified gravity behave different than in General Relativity, it 
could give a signature of the presence of higher order terms in the gravity action, as the Gauss-Bonnet invariant, 
when structure formation is studied and simulations are performed, what should be explored in the future. 
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